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Abstract 

ZnCu3(OH)6Cl2 (S = |) is a promising new candidate for an ideal Kagome Heisenberg anti- 
ferromagnet, because there is no magnetic phase transition down to ~50 mK. We investigated its 
local magnetic and lattice environments with NMR techniques. We demonstrate that the intrinsic 
local spin susceptibility decreases toward T = 0, but that slow freezing of the lattice near ~50 K, 
presumably associated with OH bonds, contributes to a large increase of local spin susceptibility 
and its distribution. Spin dynamics near T = obey a power-law behavior in high magnetic fields. 

PACS numbers: 75.40.Gb, 76.60.-k 
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A major challenge in condensed matter physics today is identifying a model material for 
investigating spin liquid [l, 2]. Searching for exotic electronic states without magnetic long 
range order, such as Kagome Heisenberg antiferromagnets, constitutes a common thread in 
a wide range of research fields, from high temperature superconductivity to low dimensional 
quantum magnetism. Over the last decade, many candidate materials have been investi- 
g ate d a S model system for a Kagome lattice g £ fl Q. However, t h e y mostly exhibit a 
magnetically ordered or spin-glass-like state at low temperatures. A recent breakthrough in 
the hunt for a spin liquid state Q is the successful synthesis Q| and characterization Q of 
ZnCu3(OH) 6 Cl2 (herbertsmithite), a chemically pure spin S = | Kagome lattice. As shown 
in Fig. 1, three Cu 2+ ions form a triangle, and a network of corner-shared triangles form a 
Kagome lattice. The S = ~ spins on Cu sites are mutually frustrated by antiferromagnetic 
super-exchange interaction J ~ 170 K Q, 9], hence the possibility of a spin liquid ground 
state. 

Recent measurements of ZnCu3(OH) 6 Cl2 with bulk magnetic susceptibility, Yw/c js} , spe- 
cific heatjs], neutron scattering on powders 8], //SR10, and 35 C1 NMR 10] have es- 
tablished that ZnCu3(OH)eCl2 remains paramagnetic down to at least ~ 50 mK with no 
evidence of magnetic long range order. These findings indeed point towards the possible 
realization of a frustrated spin liquid state with the Kagome symmetry. However, very little 
is known beyond the paramagnetic nature of the ground state. For example, the bulk aver- 
aged susceptibility, Xbuik, reveals a mysterious sharp increase below ~ 50 K 8]. This clearly 
contradicts the predictions of various theoretical calculations: series expansions predict a 
decrease of Xbuik below T ~ J/6 with a gap fl^l, while the recent Dirac Fermion model 
predicts linear behavior in T towards T = [14]. Does this apparent contradiction mean 
that ZnCu3(OH) 6 Cl2 is not a good Kagome model system after all, or that extrinsic effects 
other than Kagome Heisenberg interaction, such as mixing of Zn (S = 0) into Cu (S = k) 
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sites|15l. Il6j and Dzyaloshinsky-Moriya (DM) interactions |9l. Il7|. simply mask the intrinsic 
Kagome behavior below ~ 50 Kl What about spin dynamics? Do spin fluctuations slow 
down toward a critical point, or are they gapped 12j, 



In this Letter, we report a 63 Cu, 35 CI, and X H NMR investigation of ZnCu 3 (OH) 6 Cl 2 for 
a broad range of magnetic fields and frequencies. Taking full advantage of the local nature 
of NMR techniques, we uncover hitherto unknown properties of ZnCu3(OH) 6 Cl2. First, 
from the observation of the broadening of 35 C1 NMR lineshapes, we will demonstrate that 
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local spin susceptibility, xioc, has a large distribution throughout the sample. Moreover, the 
smallest components of xioc actually saturate and even decrease with T below T ~ 0.2J, 
even though the bulk averaged Xbulk increases as ~ 1/T. The observed decrease of xioc is 
precisely what the intrinsic spin susceptibility of a Kagome Heisenberg antiferromagnet is 
expected to show. Second, from the comparison of 1 H and 35 C1 nuclear spin-lattice relaxation 
rates, 1,35 (1/Ti), we present unambiguous evidence for slow freezing of the lattice near ~ 
50 K, most likely due to orientational disorder of OH bonds. We suggest that this subtle 
freezing of lattice distortion enhances the DM interactions, and is key to understanding the 
aforementioned upturn of bulk-averaged spin susceptibility, Xbuik, below ~ 50 K that masks 
the intrinsic Kagome behavior of xioc- Third, we demonstrate from the measurements of the 
63 Cu nuclear spin-lattice relaxation rate, 63 (1/Ti), that in the presence of a high magnetic 
field the low frequency Cu spin fluctuations grow without a gap below ~ 30 K satisfying a 
simple power-law. 

In Fig. 2, we show representative 35 C1 NMR lineshapes. For these measurements, we 
cured a powder sample in glue in a magnetic field of 9 Tesla. From powder x-ray diffraction 
measurements, we confirmed that approximately 20% of the powder is uniaxially aligned 
along the c-axis. In fact, we observe a sharp c-axis central peak near 35.02 MHz (marked 
as B//c in Fig. 2) arising from particles oriented along the c-axis. The "double horns" 
marked as #1 and #2 are split by the nuclear quadrupole interaction, and arise from the 
randomly oriented portion of the powder (i.e. 80% of the sample) [lOj. Notice that the 
whole 35 C1 NMR lineshape begins to tail-off toward lower frequencies below ~50 K. The 
resonance frequency of the sharp c-axis central peak and its distribution depends on the 
NMR Knight shift, 35 K, induced by xioc- Hence the observed line broadening implies that 
Xioc varies depending on the location within the sample below ~ 50 K. 

In Fie;. 3, we summarize the 35 C1 NMR Knight shifts 35 if and 35 K 1/2 deduced from the 
lineshapes, together with Xbuik as observed by SQUID. 35 if corresponds to the central peak 
above ~ 45 K as determined by FFT techniques. Below ~ 45 K, where the central peak 
is smeared out by line broadening, we determined 35 if as the higher frequency edge of the 
central peak from point-by-point measurements, i.e. 35 if represents the smallest component 
of the distributed xioc- 35 -^1/2 corresponds to the half-intensity position of the central peak 
on the lower frequency side of the spectrum. Quite generally, 35 if = A h fXi OC ai + 35 K chem , 
where Ahf is the magnetic hyperfine interaction between 35 CI nuclear spin and nearby Cu 
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electron spins, and 35 K c h em is a very small, temperature independent chemical shift. In the 
present case, from the comparison with Xbulk, we can estimate Ahf ~ —3.7 ± 0.7 kOe/nB- 
The negative sign of Ahf makes the overall sign of 35 K negative. Accordingly, we have 
inverted the vertical scale of Fig. 3. 

We wish to comment on two important aspects of Fig. 3. First, 35 K follows Curie- Weiss 
behavior all the way from 295 K down to ~25 K. This clearly differs from Xbuik which begins 
to deviate from Curie- Weiss behavior below temperatures as high as ~150 K [8]. On the 
other hand, series-expansion calculations indicate that the Kagome lattice follows Curie- 



Weiss behavior down to T ~ J/6 ~ 25 K 12]. Our Knight shift data demonstrate that 



Xiocai of some ideal segments of the Cu 2+ Kagome lattice in this material indeed show such 
behavior. Theoretical models also predict that below T ~ J/6 ~ 25 K, Xbuik begins to 



decrease exponentially with a small gap 



12], or linear ly|l4j. As shown in Fig. 2, the 35 CI 



NMR lineshape begins to transfer some spectral weight to higher frequencies below 25 K. 
Recalling that Ahf is negative, this is consistent with vanishing spin susceptibility, xiocai, 
near T = for some parts of the Kagome lattice. After the initial submission of this Letter, 
Olariu et al. also observed a similar decrease of 17 NMR Knight shift and arrived at the 
same conclusion [l8|. 

Another important aspect of Fig. 3 is that 35 Ki/2 begins to deviate from the aforemen- 
tioned Curie- Weiss fit in a manner similar to Xbuik, as the 35 C1 NMR line gradually broadens 
to lower frequencies. It is important to realize that 35 (1/Ti) also begins to increase in the 
same temperature range below ~150 K (see Fig. 4). Below 50 K, where 35 (1/Ti) shows a 
peak, the 35 C1 NMR line shows a dramatic broadening to lower frequencies. Fig. 2 and Fig. 
3 establish that 35 i^i/2 follows the same trend as Xbuiki i- e - some segments of the Kagome 
lattice have large and distributed local spin susceptibility Xiocai-, an d their temperature de- 
pendence is different from the smaller Xiocai as represented by 35 K. Xbuik simply represents 
a bulk average of xiocai- 

In passing, we recall that earlier /iSR Knight shift K^sr measurements by Ofer et al. [lOj] 
showed identical behavior between K^sr and Xbuik- They concluded that the upturn of Xbuik 
below 50 K is not caused by impurity spins but is a bulk phenomenon. Our new results 
in Fig. 3 do not contradict these /iSR data. K^sr was deduced by assuming a Gaussian 
distribution of Xiocah hence by default K^ SR represents the central value of the presumed 
Gaussian distribution. That explains why K^sr shows behavior similar to Xbuik and 35 Ki/2- 
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Next, we turn our attention to the dynamics of lattice and spin degrees of freedom. Fig. 4 
shows the temperature dependence of the 35 C1 nuclear spin-lattice relaxation rate, 35 (1/Ti), 
measured at the central peak frequency in various magnetic fields, B. We also plot 1 (1/T 1 ) 
for X H NMR in 0.9 Tesla, and 63 (1/Ti) for 63 Cu NMR in 8 Tesla. We have overlaid 1,63 (l/Ti) 
on 35 (1/Ti) measured in comparable magnetic fields by scaling the vertical axis. 

We can draw a number of conclusions from Fig. 4. First, let us focus on T and B 
independent results of 35 (1/Ti) above ~150 K. This high temperature regime is easily un- 
derstandable within Moriya's theory for the exchange narrowing limit of Heisenberg antifer- 



romagnets, where we should expect 35 (1/Ti) 

exc rsJ 

A 2 hf /J = constant for T > J ~ 170 jisj. 



l hf/ 

If we assume J ~ 170K and A h f ~ —4 kOe//i B , we can estimate 35 (1/Ti) exc ~ 4 sec -1 . This 
is in excellent agreement with our result. 

Another important feature is that 35 (1/Ti) begins to increase below ~150 K and peaks 
near ~ 50 K. Since 35 CI is a quadrupolar nucleus with nuclear spin I = |, the observed 
enhancement may be caused by slow fluctuations of the lattice via nuclear quadrupole inter- 
actions, as well as by Cu spin fluctuations. To discern the two possibilities, we show 1 (l/Ti) 
of 1 H measured at 0.9 Tesla for comparison. The 1 (l/Ti) data nicely interpolate 35 (1/Ti) in 
the field-independent regime above ~150 K and 35 (1/Ti) measured at a comparable mag- 
netic field (1 Tesla) at low-temperatures, without a peak near ~ 50 K. Since 1 H has I — |, 
1 (1/T 1 ) has no contributions from lattice fluctuations. Therefore we conclude that the peak 
o/ 35 (l/Ti) near ~ 50 K arises from enhancement of lattice fluctuations at the NMR fre- 
quency. The peak of 35 (1/Ti) shifts to progressively lower temperatures as we lower the 35 CI 
NMR frequency, from 50 K (35 MHz at 8.3 Tesla), 46 K (18 MHz at 4.4 Tesla) to 40 K 
(10 MHz at 2.4 Tesla). This means that the typical frequency scale of lattice fluctuations is 
35 MHz at 50 K, 18 MHz at 46 K, and 10 MHz at 40 K, and the Kagome lattice becomes 
static below 40 K. Given that no structural phase transition has been detected by x-ray 
and neutron scattering techniques, the observed freezing of the lattice near ~ 50 K must 
be a very subtle effect. In fact, our careful measurements of the quadrupole split ±| to ±| 
satellite transitions of 35 CI NMR didn't detect any noticeable changes, either. All pieces 
put together, we suggest that the lightest elements in the lattice, i.e. OH bonds, must be 
freezing with random orientations, with only subtle effects on heavier atoms. 

Regardless of the exact nature of the freezing of the lattice, our observation provides a 
major clue to understanding the mysterious behaviors of Xbuik and Xiocai- Recent structural 
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studies revealed that up to 6 % of Cu sites may be occupied by Zn to create unpaired defect 
spins [3]. There is no doubt that such antisite disorder would enhance Xbuik and contribute 
to the large distribution of Xiocai at low temperatures. However, it is important to realize 
that Xbuik begins to deviate from high temperature Curie- Weiss behavior below ~150 K 

8|, exactly where 35 (1/Ti) begins to grow due to slowing of lattice fluctuations toward 
~ 50 K. Furthermore, the crossover between the two different Curie- Weiss behaviors of 
Xbuik and 35 K 1 / 2 takes place precisely when the lattice fluctuations die out below ~50 K 
Evidently, the anomaly of the lattice correlates well with that of the spins, suggesting it 
must be playing a major role in the deviation of spin susceptibility from the theoretically 
expected behavior of ideal Kagome Heisenberg antiferromagnets. In fact, recent numerical 
simulations by Rigol and Singh showed that DM interactions can enhance spin susceptibility 
below ~50 K 0,0]. Our experimental finding naturally fits with their theoretical picture: 
when OH bonds freeze with random orientation, the hexagonal symmetry of the Kagome 
lattice slightly breaks down locally; this would progressively enhance the DM interaction 
from ~150 K to ~50 K, hence leading to an enhancement of Xbuik and 35 i^i/ 2 , as well as a 
large distribution of Xiocai- 

Finally, let us focus on the low temperature regime below 20 K, where only spin fluc- 
tuations contribute to (1/Ti). 35 (1/Ti) at ~4 K is independent of magnetic field from 
B ~0.9 Tesla up to B ~2 Tesla, but larger magnetic fields suppress 35 (1/Ti). This sug- 
gests that unpaired paramagnetic spins are fluctuating in low magnetic fields, but applied 
magnetic field decouples the weak coupling between them. The exact origin of these para- 
magnetic spins is not clear, but one obvious possibility is unpaired free spins induced in 

;he near neighbor sites of Zn (S = 0) ions occupying the Cu sites, i.e. antisite disorder 
In fact, after the initial submission of the present work, Olariu et al. reported that 
approximately 30 % of the integrated intensity of 17 NMR signals is split off from the main 
line, which they attributed to the contribution of the 4 nearest neighbor O sites of Zn ions 



occupying the Cu sites 



181 ] . Our observation of field dependence is also consistent with a 



recent report that the application of B ~ 8 Tesla suppresses defect contributions to dc spin 



susceptibility in the low T regime 



id- 



All the 63 ' 35 ' 1 (1/T 1 ) data measured at various magnetic fields and frequencies decrease 
towards the zero temperature limit. That is, we observe no hint of critical slowing down of 
spin fluctuations towards T ~ 0, hence any magnetic critical point, if it exists, is still very far 
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from our temperature reg ion. This is in agreement with an earlier report based on a limited 



set of 35 C1 NMR data 



lOj , but in remarkable contrast with earlier NMR works on frustrated 



spin systems with Kagome or triangular lattice geometry, where (1/Ti) data always show 
evidence for a magnetic long range order or freezing 6|. At B ~ 8 Tesla, both 35 (1/Xi) and 
63 (1/Ti) decrease toward T = 0, obeying a simple power law, 63 (1/Ti) = T v with rj ~ 0.5. 
Since the slope of the log-log plot of increases somewhat at higher fields, r] ~ 0.5 may 

be somewhat underestimated. We note that spin fluctuation susceptibility is proportional to 
63 (l/TiT) = T" ? ~ 1 , hence spin fluctuations grow toward T = 0. This is inconsistent with the 
exponential decrease expected for a gapped Kagome lattice in some theoretical scenarios. 
On the other hand, recent field theoretical calculations based on the Dirac Fermion model 
predicted a power law, 63 (1/Ti) = T v , with unknown critical exponent 77 |l4J. We note 
that the issue of power law behavior of spin dynamics is at the focus of the field theoretical 
approach toward spin liquid systems. It remains to be seen if the Dirac Fermion model 
would be consistent with rj > 0.5. 

To summarize, we have presented a site-by-site picture of the new Kagome material 
ZnCu3(OH)eCl2 using NMR techniques. Our NMR data revealed that both local spin sus- 
ceptibility and spin dynamics (in high magnetic fields) show aspects that are consistent 
with theoretical conjectures for ideal Kagome antiferromagnets, despite perturbations from 
lattice freezing and paramagnetic defects. 
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communications. The work at McMaster was supported by NSERC and CIFAR. We thank 
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BMB with a Graduate Fellowship Award. 
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Fig.l 

Left : A Kagome lattice. Right : Cu 2+ Kagome layer in ZnCu3(OH) 6 Cl 2 . Cl~ site is above 
or below the middle of the Cu 2+ triangles. 

Fig.2 

35 CI NMR lineshapes of the I z = +| to — \ central transition in 8.4 Tesla in a partially 
(~20 %) uniaxially aligned powder sample. The sharp peak near 35.02 MHz marked as 
"5//c" originates from the particles whose c-axis is aligned along the external magnetic 
field. Vertical lines specify the c-axis peak and edge corresponding to 35 K. Vertical dashed 
lines specify the frequency corresponding to 35 Ky 2 . 

Fig.3 

35 C1 NMR Knight shift (left scale). 35 K (blue circles). 35 K 1/2 (red triangles). The blue 
solid line is a fit to Curie- Weiss behavior, 35 K = (22 ± 7)/(T - 9 C w) + 35 K chem , where 
Qcw = —240 ± 80K. The constant background 35 K chem = 0.018% is shown by a dashed line. 
Xbuik measured by SQUID (dotted line) is also overlaid (right scale). 

Fig.4 

Temperature dependence of 35 C1 NMR spin-lattice relaxation rate 35 (1/Ti) at various 
magnetic fields (filled symbols). Solid line represents a fit to a power law, 35 (1/Ti) = T v 
with 7] = 0.47 (8.3 T), 0.44 (4.4 T), 0.2 (2.4 T and 1.0 T). X H relaxation rate in low field 
(0.9 T), ^l/Ti), and 63 Cu relaxation rate in high field (8 T), 63 (1/Ti), are also superposed 
on 35 (1/Ti) measured in comparable magnetic field. 



9 



1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 

B//c 




34.6 34.7 34.8 34.9 35 35.1 35.2 

F [MHz] 



-0.50 



-0.40^ 



-0.30-V 



-0.20- 



-0.10- 



0.00 



-i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — r 



SQUID 




0.007 



cr 



CD 



50 100 150 200 250 300 

T[K] 



-i 1 1 1 — i — i — r~r 



~i 1 1 1 — i — i — i — r 



• 8.3 T (34.6MHz) 
■ 4.4T (18.4MHz) 
▼ 2.4T (10MHz) 

♦ 1.5T (6.3MHz) 
a LOT (4.2MHz) 




10' 



OJ 



(8T, 93MHz) 
a 1-H (0.9T, 38MHz) 



* 



CO 
0) 

o 



: 1000 



10 



100 



T[K] 



